Atom trapping and guiding with a subwavelength-diameter optical fiber 
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We suggest using an evanescent wave around a thin fiber to trap atoms. We show that the gradient 
force of a red-detuned evanescent- wave field in the fundamental mode of a silica fiber can balance the 
centrifugal force when the fiber diameter is about two times smaller than the wavelength of the light 
and the component of the angular momentum of the atoms along the fiber axis is in an appropriate 
range. As an example, the system should be realizable for Cesium atoms at a temperature of less 
than 0.29 mK using a silica fiber with a radius of 0.2 /im and a 1.3-/im- wavelength light with a 
power of about 27 mW. 

PACS numbers: 32.80.Pj,32.80.Lg,03.65.Ge 



A current-carrying metallic wire and a photon-carrying 
dielectric wire are the simplest configurations for trap- 
ping and guiding cold particles. Blumel and Dietrich [j| 
considered the possibility of binding very cold neutrons 
through a magnetic trapping potential created by a thin 
wire with a current. Similar waveguides have been pro- 
posed and lately realized for binding cold neutral atoms 
0, El ■ The optical force from a dielectric wire carrying 
photons can also be used to trap and guide cold atoms 
Such an atom guide is made from a hollow opti- 
cal fiber with light propagating in the glass and tuned 
far to blue of atomic resonance pi, El- Inside the fiber, 
the evanescent wave decays exponentially away from the 
wall, producing a repulsive potential which guides atoms 
along the center axis. Alternatively, a red-detuned light 
in the hollow center of the fiber can also be used to guide 
atoms 0. In several experiments 0, cold atoms have 
been trapped and guided inside a hollow fiber where the 
optical dipole force confined them on axis. 
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FIG. 1: Schematic of atom trapping and guiding around an 
optical fiber. 



In this Rapid Communication, we present a method for 
trapping and guiding neutral atoms around a thin optical 
fiber. Our scheme is based on the use of a subwavelength- 
diameter silica fiber with a red-detuned light launched 
into it. The light wave decays away from the fiber wall 
and produces an attractive potential for neutral atoms. 
The atom trapping and guiding occur in the outside of 
the fiber. To sustain a stable trapping and guiding, the 
atoms have to be kept away from the fiber wall. This 
can be achieved by the centrifugal potential barrier. The 
centrifugal potential barrier can compensate all the po- 
tentials which diverge less rapidly than r~ 2 as r — > 0. 
We show that this can be achieved only when the fiber 
diameter is sub wavelength. (Nowadays thin fibers can be 
produced with diameters down to 50 nm Q.) 

Note that one can manipulate a very thin silica fiber 
using taper fiber technology. The essence of the tech- 
nology is to heat and pull a single-mode optical fiber 
to a very thin thickness maintaining the taper condition 
to keep adiabatically the single-mode condition [9]. The 
thin fiber has an almost vanishing core. Hence, the re- 
fractive indices that determine the fiber modes are the 
cladding refractive index nf and the refractive index no 
of the medium surrounding the fiber. The cladding index 
nf will be referred as the fiber refractive index hence- 
forth. 

Consider an atom moving around an optical fiber, see 
Fig. ^ We assume that the potential U of the atom is 
cylindrically symmetric, that is, U depends on the ra- 
dial distance r from the atom to the fiber axis z, but not 
on two other cylindrical coordinates if and z. Due to 
this symmetry, the component L z of the angular mo- 
mentum of the atom is conserved. In the eigenstate 
problem, we have L z — hm, where m is an integer, 
called the rotational quantum number. The centrifu- 
gal potential of the atom is repulsive and is given by 
U ci = h 2 (m 2 - l/4)/2Afr 2 . The radial motion of the 
atom can be treated as the one-dimensional motion of a 
particle in the effective potential U e s = U c f + U . 

There exist stable bound states for the atom if U c ^ has 
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a local minimum at a distance r = r m outside the fiber. 
This may happen only if the potential U is attractive, 
opposite to the centrifugal potential U c f. To produce such 
a potential U, we send an optical field through the fiber. 
This field generates an evanescent wave around the fiber, 
whose steep variation in the transverse plane leads to a 
gradient force on the atom. We assume that the atom is 
initially in the ground state and the detuning A of the 
field from the dominant atomic line is large compared to 
the Rabi frequency f2 and the linewidth 7. Then, the 
optical potential of the gradient force is given by [l(J, 113 
U = TiVt 2 / A. We choose a red detuning (A < 0) for the 
field to make U an attractive potential. 

We assume that the fiber is sufficiently thin that it 
has a vanishing core and it can support only a single, 
fundamental mode HEn [l^. The fiber mode charac- 
teristics are then determined by the fiber radius a, the 
light wavelength A, the fiber refractive index rif, and the 
refractive index n of the surrounding environment. In 
the linear-polarization approximation, the spatial depen- 
dence of the amplitude of the field outside the fiber is 
described by the modified Bessel function Ko(qr). Here 
the parameter q = 1/A is the inverse of the character- 
istic decay length A of the evanescent-wave field and is 
determined by the fiber eigenvalue equation [E|. Then, 
the optical potential outside the fiber can be written as 
U = -GK$(qr), where G = -fiQ 2 / AK^qa) is the cou- 
pling constant for the interaction between the evanescent 
wave and the atom. Here fl a is the Rabi frequency of the 
field at the fiber surface. 

It should be noted here that the field distribution E(r) 
corresponding to the fundamental mode HEn of the fiber 
has three nonzero components E r , E v , and E z , which 
have azimuthal variation and therefore are not cylindri- 
cally symmetric. The cylindrical symmetry of the field in 
the fundamental mode appears only in the framework of 
the linear-polarization approximation for weakly guided 
modes. Although this approximation is good for conven- 
tional fibers, it may be questionable for thin fibers since 
nf and no are very different from each other. A simple 
general way to produce a cylindrically symmetric optical 
potential is to use a circularly polarized light. The time 
average of the potential of such a field is cylindrically 
symmetric on the slow time scale of the atomic center- 
of-mass motion. 

The effective potential for the radial motion of the 
atom in the optical potential U can be written in the 
form 
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where 6 ICC — (hk) 2 /2M is the recoil energy and g = 
G/8 le c is the normalized coupling parameter. Here M is 
the mass of the atom and k is the wave number of the 
field. At the local minimum point r m , the derivative of 
U e s(r) is zero. Hence, we find r m — x m /q, where 
a solution of the equation 



with f(x) = x 3 K (x)K 1 (x) and Ad = (to 2 - l/A)q 2 /gk 2 . 
The function f(x) achieves its peak value f c = 0.2545 
at x c = 0.9331. The condition for the existence of r m is 
Ad < f c - When this condition is satisfied, Eq. (J2J has 
two solutions. The smaller one corresponds to the local 
minimum point r m of U e g(r) and the larger one corre- 
sponds to a local maximum point. It follows from the 
relation x m < x c that r m < x c A = 0.9331 A. Thus the 
distance from the local minimum point r m of the effective 
optical potential to the fiber axis is always shorter than 
the decay length A of the evanescent wave. The com- 
parison of the requirement r m > a with the inequality 
r m < x c A yields 



qa < x c £3 0.9331 



(3) 



or, equivalently, a < x c K = 0.9331 A. As shown below, 
the condition (0 can be satisfied only when the size pa- 
rameter ka of the fiber is small enough. On the other 
hand, the requirement r m > a can be satisfied only when 
f a < Ad, where f a — f(x a ) with x a — qa. Combining the 
conditions f a < Ad and Ad < f c and using the explicit 
expression for Ad, we obtain 



/„ < (to 2 - 1/4)^ < f c . 



(4) 



Thus stable bound states of the atom may exist only 
if TO min < to < TO max , where m min and TO max are 
the integer numbers closest to \f f a gk 2 / q 2 + 1/4 and 
V fcgk 2 /q 2 + 1/4, respectively. Note that an increase in 
to leads to an increase in the position r m as well as a 
decrease in the depth —U c g(r m ) of the local minimum 
of the effective optical potential. In addition, the depth 
—U c s(r m ) increases with increasing the coupling param- 
eter g. 

We calculate the decay parameter q for the evanes- 
cent wave of the fundamental mode HEn by solving the 
exact eigenvalue equation for the fiber [E|. We show 
in Fig. EJa) the characteristic dimensionless parameters 
A/ A and qa as functions of the fiber size parameter a/ A. 
We find from the figure that the condition Q is satisfied 
when 



a/A < 0.283, 



(5) 



f(x) = M, 



(2) 



that is, when the fiber is thin compared to the wavelength 
of the trapping light. The penetration length A is plotted 
in Fig. [3(b) as a function of the light wavelength A for 
various values of the fiber radius a. In these calculations, 
we took into account the dispersion of the silica-glass 
fiber index n/. The figure shows that the condition © 
is satisfied when A > 0.72, 1.06, 1.41, 1.75, and 2.09 /xm 
for a — 0.2, 0.3, 0.4, 0.5, and 0.6 /xm, respectively. 

When the atom is near to the fiber surface, we must 
take into account the van der Waals force. The van der 
Waals potential of an atom near the surface of a cylindri- 
cal dielectric rod has been calculated by Boustimi et al. 
[l3| . We use their general formula to calculate the van 
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FIG. 2: (a) Normalized penetration depth A/A and normal- 
ized decay rate qa of the evanescent wave as functions of the 
fiber size parameter a/A. The parameters used: n/ = 1.45 
and no = 1. (b) Penetration length A of the evanescent wave 
against the wavelength A of the light field for various values of 
the fiber radius a. The dispersion of the silica-glass fiber in- 
dex n/ is taken into account, and the environment refractive 
index is no = 1. 
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FIG. 3: van der Waals potential V of a ground-state cesium 
atom in the outside of a thin cylindrical silica fiber. 



der Waals potential of a ground-state cesium atom near 
a cylindrical silica fiber. We plot the results in Fig. [3] 
As seen, when r is fixed, a larger fiber radius a leads to 
a deeper and steeper van der Waals potential. When the 
van der Waals force is substantial, the effect of the cen- 
trifugal potential may be reduced and, consequently, the 
local minimum point r m of the effective optical potential 
U c s may be washed out in the total effective potential 
Utot — U c tf + V. To reduce the effect of the van der 
Waals force, we need a thin fiber. Around such a fiber, 
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FIG. 4: Total effective potential Utot of a ground-state cesium 
atom in the outside of a silica fiber with the radius of 0.2 /im. 
The wavelength and the coupling parameter of the trapping 
light are A = 1.3 /mi and g = 5330, respectively. The recoil 
energy is 6> roc = 888 Hz = 42.6 nK. 



the total effective potential E/tot may have a local mini- 
mum point near to r m . On the other hand, we can also 
reduce the effect of the van der Waals force by increasing 
r m . For a fixed value of the coupling parameter g, we 
can increase r m by considering orbits with larger values 
of the rotational quantum number m. 

We now perform numerical calculations for the total 
effective potential Utot of a ground-state cesium atom 
in the outside of a thin cylindrical silica fiber. To trap 
atoms, we use light with the wavelength A = 1.3 /im. 
The fiber radius is chosen to be a = 0.2 /im = 0.154 A, 
which is small enough to satisfy the condition J5j and 
to reduce the effect of the van der Waals force. The 
penetration length of the trapping light from the fiber 
is found to be A = 2.42 /im = 12.1a. The coupling 
parameter g should be large to produce a strong dipole 
force leading to an effective optical potential with a deep 
local minimum point. For the calculations, we choose 
g = 5330. Then, the condition J3J yields the lower and 
upper boundary values m nim — 114 and m max = 430 
for the range of m where the effective optical U e g has 
a local minimum point r m outside the fiber. For the 
calculations, we choose m = 220, 230, and 240. These 
values are not only in the interval (m m i n , m max ) but also 
are sufficiently large that r m is far away from the range 
of substantial action of the van der Waals force. We then 
expect that the total effective potential Utot has a local 
minimum point near to the point r m . 

We plot in Fig. 0] the results of our calculations for 
Utot- As seen, the total effective potential has a deep 
minimum point, which is located more than 0.4 /im far 
away from the fiber axis, not only well outside the fiber 
but also outside the range of substantial action of the 
van der Waals force. We observe that an increase in the 
rotational quantum number m leads to an increase in the 
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FIG. 5: Bound states u n for the first five levels (n = 0, 1, 
2, 3, and 4) of the radial motion a cesium atom in the total 
effective potential C/tot- The rotational quantum number is 
m = 230. All the other parameters are the same as for Fig. 

H 

position r' m of the local minimum point and a decrease in 
its depth — Ut t{r' m ). We note that, in the region of r ~ a, 
the shape of C/tot is similar to that of the attractive van 
der Waals potential V. However, in the region of r 3> a, 
where V is weak, C/tot practically coincides with C/ e ff. 

The existence of a deep local minimum of the potential 
C/tot leads to the existence of bound states of the atom. 
In Fig. [5J we plot the wave functions u n for the first five 
energy levels of the radial motion of a cesium atom in 
the potential C/tot- Since the potential is deep, the lower 
levels practically do not depend on the van der Waals 
potential. They are mainly determined by the effective 
optical potential C/ e g. However, the upper levels, which 
are not shown in the figure, are sensitive to the van der 
Waals potential and also to the boundary condition at the 
fiber surface. Tunneling from such highly excited levels 
into the narrow potential well between the repulsive hard- 



core potential and the attractive van der Waals potential 
may occur. 

The experimental realization of bound states requires 
the binding energy of atoms to be larger than their typ- 
ical thermal kinetic energy. In the case of Fig. |SJ the 
energy of the ground state is E = — 6716 # rcc . Since 
Utot{f — °°) = 0, the binding energy of an atom in the 
ground state is Eb = C/tot (f = oo) — Eq — 6716 # roc = 
0.29 mK. Thus we can trap cesium atoms around the 
fiber at a temperature of less than 0.29 mK. 

The tangential component of transverse velocity of an 
atom spinning around the fiber at the distance r m is given 
by vt — L z /Mr m = hm/Mr m . For a cesium atom spin- 
ning in a bound state with m — 230 and r m = 0.5 ^tm, we 
find vt — 0.22 m/s. The kinetic energy of this rotational 
motion is Mv T /2 = 1 mK. This energy is larger than the 
binding energy and consequently than the allowed value 
of the average thermal kinetic energy estimated above. 

We note that the coupling parameter g — 5330 corre- 
sponds to the field intensity of 1 MW/cm 2 at the fiber 
surface with the radius of 0.2 /im. Such an intensity is 
large, but it can be easily generated since the fiber is 
thin. In fact, the power of light corresponding to such an 
intensity is estimated to be about 27 mW, a very ordi- 
nary value for the laser light with the wavelength of 1.3 
fim. 

In conclusions, we have suggested and analyzed a 
method for trapping atoms based on the use of an evanes- 
cent wave around a thin silica fiber. We have shown that 
the gradient force of a red-detuned evanescent-wave field 
in the fundamental mode of a silica fiber can balance the 
centrifugal force when the fiber diameter is about two 
times smaller than the light wavelength and the compo- 
nent of the angular momentum of the atoms along the 
fiber axis is in an appropriate range. Since optical fibers 
offer the possibility of engineering evanescent fields, the 
system can be used to store, move, and manipulate cold 
atoms in a controlled manner. The system may prove 
useful as a waveguide for matter waves. 
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